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Abandoning the perfect fluid hypothesis, we investigate here the possibility that the dark energy 
equation of state (EoS) w is a nonlinear function of the energy density p. To this aim, we consider 
four different EoS describing classical fluids near thermodynamical critical points and discuss the 
main features of cosmological models made out of dust matter and a dark energy term with the 
given EoS. Each model is tested against the data on the dimensionless coordinate distance to Type 
la Supernovae and radio galaxies, the shift and the acoustic peak parameters and the positions 
of the first three peaks in the anisotropy spectrum of the comic microwave background radation. 
We propose a possible interpretation of each model in the framework of scalar field quintessence 
determining the shape of the self interaction potential V[cj>) that gives rise to each one of the 
considered thermodynamical EoS. As a general result, we demonstrate that replacing the perfect 
fiuid EoS with more generar expressions gives both the possibility of successfully solving the problem 
of cosmic acceleration escaping the resort to phantom models. 

PACS numbers: 98.80.-k, 98.80.Es, 97.60.Bw, 98.70.Dk 



I. INTRODUCTION 

The end of the 21st century has left as unexpected 
legacy a new picture of the universe depicted as a spa- 
tially flat manifold with a subcritical matter content 
presently undergoing a phase of accelerated expansion. 
An impressive amount of astrophysical evidences on dif- 
ferent scales, from the anisotropy spectrum of cosmic mi- 
crowave background radiation (hereafter CMBR) QQiQ 
to the Type la Supernovae (hereafter SNela) Hubble di- 
agram 0, Q , the large scale structure l^f and the mat- 
ter power spectrum determined by the Lya forest data 
, represent observational cornerstones that put on firm 
grounds the picture of the universe described above. 

Although the classical cosmological constant A ^ rep- 
resents the best fit to the full set of observational data 

llo| ) the well known coincidence and fine tuning prob- 
lems have lead cosmologists to look for alternative can- 
didates that are collectively referred to as dark energy. 
In the most investigated scenario, dark energy originates 
from a scalar field dubbed quintessence, running down 
its self interaction potential V{(f) so that an effective fluid 
with negative pressure contributes to the energy budget 
of the universe (for comprehensive reviews see, for in- 
stance, ^1). It is, however, also possible that dark en- 
ergy and dark matter are actually two different manifes- 
tation of the same substance. In such models, collectively 
referred to as unified dark energy (UDE), a single fluid 
with an exotic equation of state behaves as dark energy 
at the lowest ener gy s cales and as and dark matter at 
higher energies p^ll3l IT^ . It is worth remembering that 
a variant of UDE models has recently been investigated 
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by introducing models which are able to give rise to both 
inflation and cosmic acceleration llSlflal also solving the 
problem of phantom quintessence [TJ ■ 

Notwithstanding the strong efforts made to solve this 
puzzle, none of the proposed explanations is fully satis- 
factory and free of problems. This disturbing situation 
has motivated much interest toward a radically differ- 
ent approach to the problem of cosmic acceleration. It 
has therefore been suggested that cosmic speed up is an 
evidence for the need of new physics rather than a new 
ffuid. Much interest has then been devoted to models 
according to which standard matter is the only physical 
ingredient, while the Friedmann equations are modified, 
possibly as a consequence of braneworld scenarios 1^. 
In this same framework, higher order theories of gravity 
represent a valid alternative to the dark energy approach. 
Also referred to as curvature quintessence, in these mod- 
els, the gravity Lagrangian is generalized by replacing the 
Ricci scalar curvature R with a generic function f{R) so 
that an effective dark energy- like fluid appears in the 
Friedmann equations and drives the accelerated expan- 
sion. Different models of this kind have been explored 
and tested against observations considering the two pos- 
sible formulations that are obtained adopting the metric 
iil9, 20, 21, 22, 23 or the Palatini J6j formulation. 

From this overview of the different theoretical models 
proposed so far, it is clear that rather little is definitively 
known on the nature and the fundamental properties of 
the dark energy even if some model independent con- 
straints on its present day value and on its first derivative 
may be inferred from nonparametric analyses (see, e.g., 
[231 )■ It is worth noting, however, that, in all the models 
considered so far (with the remarkable exception of UDE 
models), it has been aprioristically assumed that dark 
energy behaves as a perfect fluid so that its EoS is linear 
in the energy density. Actually, from elementary ther- 
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modynamics, we know that a real fluid is never perfect 
|28| and, on the contrary, such an assumption is more 
and more inadequate as the fluid approaches its ther- 
modynamical critical points or during phase transitions. 
Given our fundamental ignorance about the properties 
of the dark sector, we cannot exclude the possibility that 
the universe is in a sort of critical point so that its con- 
stituents cannot be treated as perfect fluids. While the 
matter term, whatever its nature, may be safely mod- 
elled as a dustlike component (i.e., its EoS simplifies to 
p = 0), forcing the dark energy to be a perfect fluid is 
a rough simpliflcation that may lead to neglect the im- 
pact on the dynamics of its true properties. Moreover, 
such an unmotivated approach could lead to systemati- 
cally wrong results and hence to misleading inferences on 
the dark energy nature. Motivated by these considera- 
tions, it is therefore worth exploring what are the conse- 
quences of abandoning the perfect fluid EoS. A first step 
in this direction has been performed by Capozziello et 
al. jS^ who have considered a model in which a single 
fluid with a Van der Waals EoS accounts for both dark 
matter and dark energy (see also [s^l). From classical 
thermodynamics, we know that the Van der Waals EoS 
is best suited to describe the behaviour of real gas with 
a particular attention to the phase transitions phenom- 
ena. Actually, the Van der Waals EoS is only one of the 
possible choices in these regimes. Elaborating further on 
the idea put forward by Capozziello et al., we explore 
here other thermodynamical EoS all sharing the proper- 
ties of having been proposed to work well also for fluids in 
critical conditions. Moreover, these EoS contain the per- 
fect fluid EoS as a limiting case thus representing useful 
and more realistic generalizations. It is worth stressing 
that such an approach better reflects our ignorance of the 
dark energy nature and should prevent us from deducing 
theoretically biased conclusions on its nature. 



The plan of the paper is as follows. In Sect. II, we 
present the EoS considered giving their expressions and 
characterizing parameters. The dynamics of cosmolog- 
ical models comprising dust matter and a dark energy 
fluid with such an EoS is discussed in Sect. Ill where 
we determine the redshift evolution of the main physical 
quantities of interest. Matching with the data allows to 
investigate the viability of the different EoS and constrain 
their parameters. The method we use and the results we 
get are presented in Sect. IV, while the position of the 
peaks in the CMBR anisotropy spectrum is evaluated in 
Sect. V and compared with the WMAP determination. 
In Sect. VI, we reinterpret the models proposed in the 
framework of scalar field quintessence determining the 
self interaction potential that gives rise to a dark energy 
model with the given EoS. Conclusions are presented in 
Sect. VII, while in Appendix A we give some further de- 
tails on the EoS from the thermodynamic point of view. 



II. EQUATIONS OF STATE 

The dynamical system describing a Friedmann- 
Robertson- Walker (FRW) cosmology is given by the 
Friedmann equations ||3l| : 
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(pa/ + Px + 3px) 



SttG 



[PM + Px) , 



(1) 



(2) 



and the continuity equations for each of the two fluids : 



3i7 (a+pO = 



(3) 



where a is the scale factor, H = a/a the Hubble parame- 
ter, the dot denotes the derivative with respect to cosmic 
time and we have assumed a spatially flat universe in 
agreement with what is inferred from CMBR anisotropy 
spectrum 

HHH' Eqs.lP), IP and © are derived by 
the Einstein fleld equations and the contracted Bianchi 
identities assuming that the source of the gravitational 
fleld is a a mixture of matter with energy density pM 
and pressure pm = and an additional negative pressure 
fluid (which is usually referred to as dark energy) with 
energy density px and pressure px ■ To close the system 
and determine the evolution of the scale factor a and 
of the other quantities of interest, the equation of state 
(hereafter EoS) of the dark energy fluid (i.e. a relation 
between px and px) is needed. 

Unfortunately, this is a daunting task given our com- 
plete ignorance of the dark energy nature and of its fun- 
damental properties. Motivated by the discussion in the 
introduction, we explore here some EoS all sharing the 
properties of working well even when the fluid is near 
critical points or phase transitions. A textbook example 
is the Van der Waals EoS : 



Px 



IPX 

l-Ppx 



(4) 



wherec a and (3, in the thermodynamics analogy, may be 
related to limiting values of the pressure and the volume, 
while 7 is the usual barotropic factor. The Van der Waals 
EoS reduces to the perfect fluid case in the limit a, /3 ^ 0. 
The dynamics of the corresponding cosmological model 
has been yet investigated both in the framework of UDE 
models |23 and as a dark energy fluid 30] so that we 
do not consider it again here. On the other hand, there 
arc other EoS that are worth investigating. However, 
we limit our attention to EoS described by two param- 
eters only in order to both narrow the parameter space 
to explore and avoid introducing too large a degeneracy 
among the quantities we have to determine. 

Let us first define 77(z) and p{z) as px{z) / Pcrit and 
Px{z)/pcrit respectively, being pcnt = SHq/SttG the 
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present day critical density of the universe. The EoS 
may then be evaluated as w — px/px = p/v- For the 
different models we consider, w is a non linear function 
of rj and is given as follows. 

i. Redlich - Kwong : 



WRK = /3 X 



1 - - 2V2aT] 
1 - (f - V2)a'n 



(5) 



deceleration parameter q = —aajc? and evaluating the 
result at z = 0, we get the well known relation: 



90 



-S^xWq 



(10) 



Introducing one of Eqs.(|SJl - ^ into Ea. ljlUII with the con- 
dition 7/(z = 0) = fix, we can express (5 in terms of goj V 
and thus obtaining what follows. 

i. Redlich - Kwong : 



ii. Modified Berthelot : 
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(11) 



ii. Modified Berthelot: 



iii. Dieterici: 
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/?exp [2(1 - arj)] 
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(7) 



iii. Dieterici : 



(2go~l)(l + y) 
3(1 - ^m) 



(12) 



iv. Peng - Robinson : 
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(1 + ari)/{l — arj) + my 



(8) 



with 



1.487 and ct ~ 0.253. 



In Eqs.JSJ - (|H1), the two parameters a and /? are related 
to the critical values of density and pressure of the fluid. 
In particular, for all cases, a oc pcrit/Pc while /S (x pc/pc 
with pc and Pc the values of the energy density and pres- 
sure respectively at the critical point of the fluid^. Note 
that, for a = 0, all the EoS above reduces to w = est, 
i.e. to the perfect fluid one. The condition a = is 
achieved for an infinite critical density which means that 
the fluid has no critical points. This is indeed the case of 
the perfect fluid and is the reason why such a description 
is highly unrealistic. 

It is convenient, however, in the application to ex- 
press these two parameters in terms of more handable 
and meaningful quantities. To this aim, we first define 



y = aTj{z = 0) = a{l - ^m) ^a = y/{l - ^m) (9) 

where we have used the flatness condition J^m + ^x = 
1. Combining Eqs.lQJ and (0), using the definition of 



^ See the Appendix for the definition of critical points and the 
exact expression of a and 
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3(1 - ^m) 
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iv. Peng - Robinson : 



Cb(2go-l)(t/-l)(y2-2y-l) 
3(1 - nM)[cayil -y) + cbiy^ -2y- 1)] 



(14) 



Note that it is y rather than a to determine f3 so that 
it is this parameter that will be constrained by the fit- 
ting procedure. Moreover, go ^^nd Qm are more familiar 
quantities than (3 so that it is easier to choose intervals 
over which they can take values. 



III. REDSHIFT EVOLUTION 

For a given EoS, it is possible to determine how the 
main physical quantities (such as the energy density, the 
EoS and the Hubble parameter) evolves with the redshift 
z. To this aim, we first change variable from t to z in the 
continuity equation Q which thus is rewritten as : 



drj 3(1 + w)ri{z) 
dz 1 + z 



(15) 



with the initial condition ri{z = 0) = ^x- Note that, 
for the EoS we are considering, Ea. H15|l is a first order 
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FIG. 1: The evolution against the redshift of the dimension- 
less energy density ri{z) for the RK (upper left), MB (up- 
per right), Dt (lower left) and PR (lower right) EoS respec- 
tively. For all models, we set {qo,ilM) ~ (—0.5,0.3). Short 
dashed, solid and long dashed lines refer to different values of 
y, namely y = 0.5, 0.75, 1.0 for the RK model, y = 1.5, 2.5, 3.5 
for the MB one, y = 0.5, 0.7, 0.9 for the Dt and PR models. 



nonlinear differential equation that cannot be solved ana- 
lytically. However, a numerical integration is straightfor- 
ward provided that the parameters {qo, y, ^m) are given. 
Fig.n shows the results for the different EoS described 
in the previous section for some illustrative set of pa- 
rameters. Note that, hereon, to shorten the notation, we 
will use the acronyms RK, MB, Dt and PR referring to 
the Redlich- Kwong, Modified Berthelot, Dieterici and 
Peng - Robinson cases respectively. 

As it is clearly shown, not surprisingly, different EoS 
may lead to radically different evolutions for the energy 
density. This is particularly true comparing the upper 
panels with the lower ones. For the RK and MB cases, 
log 77 is an almost linear increasing function of z, i.e. r]{z) 
has an approximately power -law like decrease with the 
cosmic time t over a large range. As a consequence, in 
the far past, the dark energy component does not fade 
away, but still contributes to the energy budget during 
the usually matter dominated epoch. This behaviour 
may be problematic for the impact on structure forma- 
tion and nucleosynthesis. For the RK case, this problem 
may be particularly worrisome since, as can be inferred 
from Fig.^ for high z, rj ^ (1 + z)'^ with 7 larger than 3 
for some combination of the parameters {qo,y, S^m)- The 
situation is less dramatic for the AIB EoS since, although 
we still get 77 ~ (1-l-z)''' for high z, now 7 is smaller than 3 
so that the dark energy term becomes quite small during 
the matter dominated era. Note also that the evolution 
of r](z) only weakly depends on y for the M B model thus 
suggesting a serious degeneracy in this quantity. The sit- 
uation is radically different for the Dt and PR models. 
As Fig.Qshows, in these cases, 77(2) quickly approaches a 
constant value so that, in the past, the dark energy com- 
ponent does not disappear, but plays the same role as 
a cosmological constant term. Note, in particular, that 
this regime is achieved very soon for the PR EoS in which 




FIG. 2: The evolution against the redshift of the EoS pa- 
rameter ■w{z) for the RK (upper left), MB (upper right), Dt 
(lower left) and PR (lower right) EoS respectively. Model 
parameters are set as in Fig.0 



case r](z) is almost constant for quite small values of z. 
These results are reassuring since the energy density of 
the dark energy component becomes vanishingly small 
with respect to that of the matter during both the struc- 
ture formation and nucleosynthesis epochs so that we are 
quite confident that these processes are not altered by the 
use of unusual EoS. 

Having determined 77(2:), it is now straightforward to 
compute how the EoS depend on the redshift. The re- 
sult is shown in Fig.[21 where we plot w{z) for the differ- 
ent models adopting for the model parameters the same 
values as in Fig.Q] Although the behaviour of 77(2) is 
qualitatively similar in some cases, the shape of w{z) is 
radically different for the models we are considering so 
that we discuss them separately. 

First, let us consider the RK EoS. For the models in 
the upper left panel of Fig.|21 the EoS turns out to be 
an increasing function of the redshift 2 with the largest 
value of y leading to higher w at high 2. In particular, 
w may become positive for sufficiently large values of 
y. However, we have checked that this result strongly 
depends on flM- Actually, for values of f^M > 0.35, the 
EoS becomes more and more negative as 2 increases so 
that the fluid behaves as kind of superphantom. As a 
general rule, however, we stress that, for y < I, w{z ~ 
0) ~ —1 so that in the present epoch the EoS mimics 
that of the cosmological constant. 

The M B EoS is shown in the upper right panel of 
Fig.El but the main trend should be inferred directly 
from Eq.(|SJ|. Since 77(2) is an increasing function of 2, 
it is easy to understand that, whatever are the values 
of (go,y, ^^m), wlz) will vanish in the past so that a 
dust -like EoS is asymptotically achieved. Note, how- 
ever, that the convergence may be quite slow depending 
on y : the larger is y, the higher is w at a given z so 
that the quicker is the convergence toward the asymp- 
totic null value. Given this behaviour, it is tempting to 
use the MB EoS as a proposal for a UDE model. From 



5 



the point of view of the parameters, this model may be 
obtained imposing — ^ — ^t with Qf, the baryon den- 
sity parameter. However, we prefer to not fix fix and 
determine it later from matching with the data. 

Let us consider now w{z) for the Dt parametrization 
(lower left panel in Fig.[21). In contrast with the other 
cases considered, is not a monotonic function of the 
redshift, but it has rather an asymmetric bell -shaped 
behaviour. In particular, the height of the peak is larger 
for smaller values of y and its position shifts toward right 
(i.e., larger values of z) with the increasing of y. The 
most remarkable feature is, however, the asymptotic ap- 
proach toward the cosmological constant value w = ~1 
that is achieved later for smaller y. A similar behaviour 
is consistent with the result shown in Fig. ^ according to 
which ri{z) becomes constant for values of z sufficiently 
high. Comparing the two plots, we see that r]{z) starts 
being approximately constant as soon as w{z) is indis- 
tinguishable from — 1 so that everything works as for the 
cosmological constant. Even if not shown in the plot, we 
note that w{z) approaches —1 also in the limit z ^ —1, 
i.e. in the asymptotic future, so that a de Sitter like 
expansion is achieved. 

Finally, let us discuss the case of the PR EoS which is 
depicted in the lower right panel of Fig.[21 As a striking 
result, we get that w{z) starts from a value very close 
to -1 and very soon reaches w = —1 after which it does 
not change anymore. This behaviour nicely explains why 
ri{z) is approximately constant over almost the full evo- 
lutionary history. It is also worth noting that, although 
w{z) depends significantly on y, the numerical change in 
its value is too small to be detected. As a consequence, 
it is likely that matching with the data will be unable to 
efficiently constrain this parameter. 

Another interesting dynamical quantity is the deceler- 
ation parameter q. Combining the Friedmann equations, 
it is straightforward to get : 



the present age of the universe), it turns out that air) 
is almost linear over the most of the universe evolution 
and, what is more interesting, is independent of y. Ac- 
tually, this is only a result of having used the dimension- 
less time T rather than the physical time t. Since, as we 
have checked, tg depends significantly on the combina- 
tion of the model parameters {qo,y,flM), transforming 
from a(r) to a{t) introduce the expected dependence on 
y. As a general result, a(t) does not diverge in any finite 
time so that any Big Rip is avoided even if w may lie 
today in the phantom {w < —I) regime. 



IV. CONSTRAINING THE EOS 

The discussion above has shown that the dynamics of a 
cosmological model filled with dust matter and a dark en- 
ergy fluid whose EoS is one of those proposed in Sect. II 
is reasonable and not affected by any pathological be- 
haviour (provided the parameters are chosen with some 
care in the case of the RK model). It is therefore inter- 
esting to compare the models with the available obser- 
vations in order to both investigate the viability of the 
model itself and constrain its parameters. 



A. The method 

In order to constrain the EoS characterizing parame- 
ters, we maximize the following likelihood function : 



C oc exp 



xHp) 



(17) 



where p denotes the set of model parameters and the 
pseudo - merit function is defined as : 



w{z)rj{z) 



N 



2 2f]M(l 



T]{z) 



(16) 



so that, having yet evaluated both 77(2) and 'w{z), it is 
immediate to compute q{z). It turns out that, for all the 
EoS considered, the evolution of q{z) is quite similar over 
the redshift range probed by the most of the available 
data. Moreover, it is remarkable that there is almost no 
dependence at all on y for the MB, Dt and PR models, 
while a weak dependence is present in the case of the 
RK EoS. In this latter case, it is important to stress that 
the adopted value of flm plays a fundamental role with 
values of > 0.35 leading to q{z) < for all z > 0, 
i.e. these models are never decelerating. For all other 
cases, instead, q{z) changes sign so that the transition 
redshift, defined as q{zT) ~ 0, turns out to be positive in 
agreement with some recent estimates. 

As a final issue, we have also numerically evaluated 
how the scale factor depends on the cosmic time t. In- 
troducing the normalized time variable r = t/to (with to 



„th 



(2j,p) 



^obs 



7^(p) 



1.716 



0.062 



A{p) - 0.469 



0.017 



.(18) 



Let us discuss briefly the different terms entering Ea. (|18|l . 
In the first one, we consider the dimensionless coordinate 
distance y to an object at redshift z defined as : 



r(z) 



dz' 

W) 



(19) 



and related to the usual luminosity distance as = 
(1 -f z)r{z). Daly & Djorgovki 32J have compiled a sam- 
ple comprising data on y{z) for the 157 SNela in the 
Riess et al. Q Gold dataset and 20 radiogalaxies from 
[33I I , summarized in Tables 1 and 2 of 32j . As a prelim- 
inary step, they have fitted the linear Hubble law to a 
large set of low redshift (z < 0.1) SNela thus obtaining: 
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h = 0.664 ± 0.008 . 

We thus set h = 0.664 in order to be consistent with 
their work, but we have checked that varying h in the 
68% CL quoted above does not aher the main resuhs. 
Furthermore, the value we are using is consistent also 
with iJo = 72 ± 8 km s'^ Mpc'^ given by the HST Key 
project 13 based on the local distance ladder and with 
the estimates coming from the time delay in multiply 
imaged quasars and the Sunyaev- Zel'dovich effect 
in X - ray emitting clusters 36] . 

The second term in Eq. (|f 8() makes it possible to extend 
the redshift range over which y{z) is probed resorting to 
the distance to the last scattering surface. Actually, what 
can be determined from the CMBR anisotropy spectrum 
is the so called shift parameter defined as (37i . i38»| : 



R 



iy{zis) 



(20) 



where zis is the redshift of the last scattering surface 
which can be approximated as |3^.: 



zis = 1048 (1 + 0.00124:Uj-° '^^^) (1 



92^ 



(21) 



with LUi — flih^ (with i — b, M for baryons and total 
matter respectively) and (51,32) given in Ref. '3^. The 
parameter ujb is well constrained by the baryogenesis cal- 
culations contrasted to the observed abundances of pri- 
mordial elements. Using this method, Kirkman et al. 
[iof have determined : 

ujb = 0.0214 ±0.0020 . 

Neglecting the small error, we thus set Wf, = 0.0214 and 
use this value to determine zig. It is worth noting, how- 
ever, that the exact value of zis has a negligible impact 
on the results and setting zis = 1100 does not change 
none of the constraints on the other model parameters. 

Finally, the third term in the definition of takes into 
account the recent measurements of the acoustic peak in 
the large scale correlation function at 100 Mpc sep- 
aration detected by Eisenstein et al. fS*] using a sample 
of 46748 luminous r ed g alaxies (LRG) selected from the 
SDSS Main Sample |43|. Actually, rather than the posi- 
tion of acoustic peak itself, a closely related quantity is 
better constrained from these data defined as 41] : 



A 



ZLRG 



zlrg 



E{zlrg) 



V^{zlrg] 



1/3 



(22) 



with zlug — 0.35 the effective redshift of the LRG sam- 
ple. As it is clear, the A parameter depends not only on 
the dimensionless coordinate distance (and thus on the 
integrated expansion rate), but also on 17 j\/ and E{z) ex- 
plicitly which removes some of the degeneracies intrinsic 



in distance fitting methods. Therefore, it is particularly 
interesting to include A as a, further constraint on the 
model parameters using its measured value ^3] : 

A = 0.469 ±0.017 . 

Note that, although similar to the usual reduced in- 
troduced in statistics, the reduced (i.e., the ratio be- 
tween the Eind the number of degrees of freedom) is 
not forced to be 1 for the best fit model because of the 
presence of the priors on TZ and A and since the uncer- 
tainties cTi are not Gaussian distributed, but take care 
of both statistical errors and systematic uncertainties. 
With the definition H17|l of the likelihood function, the 
best fit model parameters are those that maximize 'C(p). 
However, to constrain a given parameter p^, one resorts 
to the marginalized likelihood function defined as : 

Cp,{pi) (X j dpi . . . j dpi^i J dpi+i... J dpnC{p) (23) 

that is normalized at unity at maximum. Denoting with 
Xq is the value of the x^ for the best fit model, the 
1 and 2(7 confidence regions are determined by impos- 
ing Ax^ = ^ Xo = 1 ^^'^ ~ 4 respectively. 



B. Results 

We have applied the likelihood analysis described 
above to the four EoS presented in Sect. II obtaining con- 
straints on the model parameters {qo, y, ^m)- The results 
obtained are summarized in Table I where we give the 
best fit values and 68% and 95% confidence ranges for 
each parameter. Note that the range tested for y is set on 
a case by case basis as discussed in the following. Given 
(go, 2/j ^m) for an EoS, we may also evaluate some other 
interesting physical quantities. Since the uncertainties 
on the model parameters are not Gaussian distributed, a 
naive propagation of the errors is not possible. Moreover, 
we have not an analytical expression for some quantities 
such as, e.g., the transition redshift. We thus estimate 
the 68% and 95% confidence ranges on the derived quan- 
tities by randomly generating 20000 points (go,2/, ^^m) 
using the marginalized likelihood functions of each pa- 
rameter and then deriving the likelihood function of the 
corresponding quantity. Although not statistically well 
motivated, this procedure gives a conservative estimate 
of the uncertainties which is enough for our aims. 

As a general result, we note that the constraints on 
both qo and ^Im turn out to be essentially model inde- 
pendent. Moreover, they are consistent with recent esti- 
mates obtained using different datasets and dark energy 
models with a perfect fluid EoS (constant or redshift de- 
pendent). In particular, both values are quite similar to 
those predicted for the concordance ACDM model yield- 
ing {qQ,il.M) — (—0.5,0.3). Actually, this is not much 
surprising. As discussed in Sect. Ill, the four EoS con- 
sidered mimic well the cosmological constant for small 
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Id 




y 


Am 


A 


7^ 


Zt 


to (Gyr) 


RK 
MB 
Dt 
PR 


n /iQn+0 005+0.009 

"•^■="-'-0.005-0.009 

n irc;+0.04+0.09 
U.OO_o.04-0.09 

n cc;+0.04+0.09 
"•^'-'-0.04-0.09 

n rr+0.07+0.32 
"•^"-•-0.07-0.46 


1 Qr,+0.02+0.04 
±.»U_Q Q2_o.o4 

(< 4.2 < 6.0) 
0.9 (< 0.9) 
0.9 (> 0.57 > 0.24) 


n occ+0. 001+0. 003 

U.OOO_Q Qm_Q QQ2 

r, 90+0.02+0.04 
"•^"-0.02-0.04 

n 90+0. 02+0. 04 
"•^"-0.02-0.04 

n 9Q+0.02+0.04 
"•^"-0.01-0.03 


n 474+0.014+0.030 
U.'±<'±_o.oi8-0.033 

0/171 +0.015+0.033 
^-0.016-0.033 

n 4 7f;+0.015+0.030 
"■^'"-0.012-0.025 


1 700+0022+0. 042 
1- loO_o.Q32-0.095 

1 700+0.021+0.045 
±. ( 'JO_o 034-0.065 

1 7/IO+0. 012+0. 026 
1- '^1^-0.012-0.027 


n CQO+0. 050+0. 094 
U.OyZ_o 042-0.082 

n 7nR+0 053+0.107 
U. (UU_o 062-0.132 

n 71 n+0. 052+0. 122 
U. 1 l"-0. 060-0. 113 


1/1 9Q+0.22+0.43 
^^^■'^^-0. 21-0. 41 

1/L ^9+0.23+0.49 
^^■■^^-0.24-0.49 

1/1 00+0.33+0.54 
-'-^■'^^-0.19-0.45 



TABLE I: Summary of the results of the hkehhood analysis of the models discussed in the text. The meaning of the entries 
is as follows. By writing x = , we mean that x is the maximum likelihood value of the considered quantity, while 

the 68% and 95% confidence ranges are {x ~ S- , x + 5+) and (x — S--, x + 5++) respectively. We use this notation to give our 
constraints on the model parameters (go, J/, ^m) and the derived quantities {A, TZ, zt, to)- For the RK case, we do not estimate 
{A,TZ,ZT,to) since the model may be discarded. For the MB and Dt case, we are able to give only upper limits on y, while 
only lower limits may be given on this same quantity for the PR model. See the text for discussion. 



values of z. Therefore, we do expect similar values for go 
and SIm since these quantities are both evaluated today 
when the difference among the concordance model and 
our ones may be hardly detected. 

The only parameter directly characterizing the differ- 
ent models is therefore y so that constraints on y are 
indeed strongly model dependent. Beside, the physical 
range for this quantity must be set on a case by case ba- 
sis thus obviously impacting the final estimate so that we 
discuss separately the results for each model. 

1 . Redlich - Kwong 

As a first issue, it is important to assess what is the 
range explored for the parameter y. Looking at Eq.(|SJ), 
it is clear that, in order to avoid unphysical divergences 
of the EoS, the condition 1 — \/3 — 2^/2ari{z) ^ must 
hold. Moreover, from Ea. l|ll|) . we get the further con- 
strain y ^ l/\/3- 2V2 in order /3 to not be diver- 
gent. We have checked (analytically and numerically) 
that choosing < y < 2 ensures that all the conditions 
quoted above are fulfilled whatever is the redshift z. Per- 
forming the likelihood analysis discussed above, we have 
obtained the constraints reported in the first row of Ta- 
ble!. Quite surprisingly, the range for y is very narrow. 
Exploring the likelihood contours in the 3 — D parame- 
ter space, we have found that there are indeed two local 
minima of the pseudo-x^ defined above. Our procedure 
select the absolute minimum thus selecting a very small 
region. There is, however, also a more subtle motiva- 
tion for the very stringent constraints obtained on y. Al- 
though not apparent from the upper left panel of Fig-Q] 
the energy density strongly depends on y in the region 
y > 1 so that also small deviations from the best fit value 
leads to significant departures from the best fitting curve 
thus leading to very strong constraints on the model pa- 



rameters. 

Although the best fit curve reproduces very well the 
data, the RK model may be excluded on the basis 
of physical considerations. Actually, for (go, 2/, hi 
the parameter space individuated by the constraints re- 
ported, the energy density increases with z faster than 
the matter one. As a result, the universe turns out to be 
ever accelerating (so that the estimated transition red- 
shift is negative) and never undergoes a matter domi- 
nated epoch in the past. Even if we have not performed 
a detailed calculation, it is nevertheless clear that such 
a situation leads to severe problems with both structure 
formation and nucleosynthesis. Note that the same qual- 
itative behaviour holds for the parameters taking values 
in the other local mimimum. Given these problems, we 
conclude that the RK EoS may be discarded and will not 
be considered anymore in the following. 



2. Modified Berthelot 

Looking at Eqs.© and H12() . it is clear that there are 
no physical motivations to impose an upper limit on j/, 
the only constraints thus being j/ > 0. In the limit 
y » 1, combining Eq.® with Eqs.Q and ((T^ . we get 
w ~ {l3/a)/ri = {2qo - 1)/[3(1 - so that the EoS 

does not depend on y in this regime. As a consequence, 
the constraints on y turns out to be quite weak and the 
likelihood analysis may put only upper limits. For the 
best fit values in the second row of Table I, Eq.® reduces 
to the perfect fiuid one with w = (250 ~ 1)/[3(1 — rijv/)] — 
— 1.12. This result could suggest that the likelihood anal- 
ysis argues in favor of no need of giving off the perfect 
fluid hypothesis. Actually, it is worth stressing that the 
marginalized likelihood function for y is quite flat so that 
values of y 7^ are perfectly viable. Indeed, the la confi- 
dence range extends up to y = 4.2 thus showing that the 
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MB EoS provides a good match with the data even when 
it significantly differs from the simplest model p = wp 
with w a constant. 

Let us now briefly comment on the possibility to use 
the MB EoS in the framework of UDE models. Should 
this approach be correct, the likelihood analysis should 
have returned ^Im — f^b, while such low values are ex- 
cluded at more than 3a level. As such, one could con- 
clude that the UDE approach may be rejected. Actually, 
one should still explore the possibility that Eq.lO is an 
effective EoS and formally decompose the energy den- 
sity px as sum of pdm and pde with the first and second 
term referring to dark matter and dark energy respec- 
tively. The EoS of this dark energy term is then eval- 
uated imposing wmb = WdePde/{pdm + Pde)- Imposing 
Pdm = ^dmPcriti^ + z)^, One should set rtdm = - 
using the value of Qm determined above and a model 
independent estimate of fit,. Investigating this scenario 
is outside our aim so that we do not speculate further on 
this interesting possibility. 

It is interesting to discuss with some detail the con- 
straints derived on some physical quantities coming from 
the likelihoods for the model parameters. As a first con- 
sistency check, we have estimated both the acoustic peak 
and the shift parameters A and TZ. Even if we have ex- 
plicitly introduced priors on them in the definition of the 
pseudo - in Eq- (|18|l , it is nevertheless possible that the 
likelihood procedure selects a region of the parameter 
space giving values of A and TZ in disagreement with the 
imposed priors^. Actually, it turns out that A and TZ 
agree very well (within la) with the measured ones, al- 
though the maximum likelihood values are slightly larger 
than the estimated ones. 

While qo < 0, the universe has entered the epoch of 
accelerated expansion only for z < zt, this latter be- 
ing the transition redshift previously defined and con- 
strained for the MB model as reported in Table I. The 
maximum likelihood value zt = 0.692 is more than 1.7a 
larger than the tentative model independent estimate of 
Riess et al. giving zt = 0.46 ± 0.13 M- It is worth not- 
ing, however, our value of zt is in good agreement with 
that predicted for the concordance ACDM model, being 
in this case zt = {2nA/nM)^^^ -1~ 0.671. This is not 
much surprising given that, in the region of the parame- 
ter space selected, the MB EoS mimics well that of the 
A term over the redshift range probed by the data. 

Finally, we consider the age of the universe obtaining 
to — 14.29 Gyr as maximum likelihood value and the 2a 
confidence range extending from 13.88 up to 14.72 Gyr. 
This result is in satisfactory agreement with previous 
model dependent estimates such as to = 13.24^q'45' Gyr 
from Tegmark et al. and to = 13.6±0.19 Gyr given by 



^ Qualitatively, this could be understood noting that the main con- 
tribution to comes from the 157 SNela, while A and TZ gives 
only a modest contribution unless the model is unreasonably dif- 
ferent from the best fit one. 



Seljak et al. [To|. Aging of globular clusters and nu- 
cleochronology |4J| give model independent (but affected 
by larger errors) estimates of to still in good agreement 
with our one. 



3. Dieterici 

Setting the range of y for the Dt EoS is a subtle task. 
Imposing that w never diverges and f3 does not vanish 
leads to the conditions y ^ 2 — [y/(l — VLM)\'r]{z) and 
J/ 7^ 2. There is, however, a further constraint motivated 
by numerical integrations of the continuity equation that 
turns out to become unstable for j/ > 1. In order to avoid 
this problem, we have searched for the constraints on y 
in the region < y < 1 only. It comes out that the 
marginalized likelihood is quite flat so that the full range 
is well within la from the best fit value. Although this 
is quite disturbing from the point of view of constraining 
the model, this is encouraging since it shows that aban- 
doning the perfect fluid EoS in favor of the Dt one still 
gives a good match with the data. As a final remark, 
let us note that the acoustic peak and the shift parame- 
ters A and TZ, the transition redshift and the age of 
the universe to estimated for the Dt model are in very 
good agreement with the same quantities obtained for 
the MB case so that we refer the reader to what already 
said above. 



J^. Peng - Rohinson 

Eq. I|14|l shows that there are two values of y such that 
the PR EoS reduces to the perfect fluid one, namely y = 

(so that a = and Eq.lO reduces to wpr = (3) and y — 

1 (giving wpB^ — 0). We have checked that values oiy > 1 
give rise to models having some pathological behaviours 
in the past (for instance, unphysical divergence of 'q{z) 
for high z) so that we have restricted our attention to 
the range < ?; < 1. Once again, the likelihood function 
is quite flat so that we are able only to give lower limits 
on y. It is noteworthy, however, that the best fit value 
is now y — 0.9, that is the PR EoS does not reduce to 
that of the perfect gas. Finally, it is worth noting that 
the constraints on A, TZ, zt and to obtained for the PR 
model agree very well with those estimated for the MB 
case so that we refer the reader to what already said. 



C. Degeneracy with the ACDM and QCDM models 

The results discussed above demonstrates that the 
MB, Dt and PR EoS give rise to cosmological models 
that are in good agreement with the considered dataset. 
On the other hand, also models with a dark energy hav- 
ing a perfect fluid EoS provides a very good match to 
the same dataset. This consideration suggests that a 
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Input Model 


ACDM 


QCDM 


T J 

id 








Wx 


MB 


1.0,0.28 


r, 9Q+0.02+0.04 
"•^°-0.02-0.04 


U.OU r\ r\A (-1 r\Q 

— U . U^: — U . UO 


1 0^+0.12+0.22 
— u.±u — u.oo 


Dt 


0.25,0.28 


r, 9f;+0.02+0.04 
"•^"-0.02-0.04 


r, n^+0.03 + 0.06 
'-'■'^^-0.04-0.08 


1 97+0.15+0.27 
^■^ ' -0.17-0.38 


PR 


0.35,0.28 


r, 97+0.03+0.06 
-0.02-0.05 


97+0.04+0.08 
^■^ ' -0.04-0.08 


Q7+0.09+0.18 
_o.i2-0.26 



TABLE II: Summary of the results of the likelihood anal- 
ysis on the simulated dataset described in the text. The 
first two columns identifies the input model giving the EoS 
id and model parameters. In particular, for all the EoS, we 
have set qo = —0.55. The third column refers to the ACDM 
model, while fourth and fifth columns are for the QCDM 
model. Maximum likelihood values and confidence ranges are 
reported using the same scheme as in Table I. 



sort of degeneracy among the different EoS should ex- 
ist. Investigating in detail this issue needs a detailed set 
of simulations in order to understand under which con- 
ditions such a degeneracy may be broken. Although this 
is outside our aim, we nevertheless provide a preliminary 
analysis that is sufficient to get an interesting feeling of 
the problem. To this aim, we implement a quite simple 
procedure. First, we select an EoS and set its character- 
izing parameters. Then, we generate a sample of SNela 
according to the theoretical luminosity distance for the 
model with the EoS chosen before. Note that the sample 
comprises the same number of SNela of the Riess et al. 
Gold sample and have the same redshift and distance 
modulus error distribution. Finally, we fit to this dataset 
the ACDM [px = -px) and the QCDM [px = wxPx) 
model. In order to render our analysis as similar as pos- 
sible to that in Riess et al. , in the second case we impose 
the prior VIm — 0.27 ± 0.04 as done in ^^J. The parame- 
ters used in the simulations and the constraints obtained 
on the A and QCDM model parameters are summarized 
in Table II. 

Some interesting lessons may be learned from this sim- 
ple exercise. First, we note that the ACDM model fits 
well in all the cases considered. Moreover, the estimated 
VIm is quite close to the input value so that no system- 
atic errors is induced on this parameter. Note that this 
result does not depend on the particular choice of the 
EoS parameters provided they lie in the confidence ranges 
summarized in Table I. Actually, such a result could be 
expected considering that, over the redshift range probed 
by the SNela sample and for the values of parameters cho- 
sen, the three EoS chosen mimic quite well the ACDM 
model so that it is not surprising that the simulated data 
may be fitted by the concordance scenario. It is still 
more interesting to consider the results from fitting the 
QCDM model to the simulated dataset. While the esti- 
mated fiji/ is again quite similar to the input value (al- 



though biased and formally not in agreement within the 
errors), the constraints on wx may extend in the phan- 
tom region {wx < —1) depending on the EoS adopted 
and the parameters chosen. This preliminary test sug- 
gests a possible way to escape the need of the problem- 
atic phantom dark energy (i.e. a negative pressure fluid 
with Wx < — !)■ Indeed, Tablell shows that wx < — 1 
may be the consequence of forcing the perfect fluid EoS 
to flt a cosmological model where the true dark energy 
EoS is not the perfect fluid one. This intriguing scenario 
have, however, to be further investigated with a more 
careful and extensive set of simulated dataset also taking 
into account other possible probes such as the priors on 
A and TZ or the gas mass fraction in galaxy clusters. 



V. THE CMBR PEAKS POSITION 

The analysis presented above has convincingly shown 
that dark energy models with EoS given by the MB, 
Dt and PR parametrizations are indeed viable alterna- 
tives with respect to the usual perfect fluid assumption. 
Indeed, the r{z) diagram, the acoustic peak A and the 
shift parameter TZ are correctly predicted and also the 
estimated age of the universe is in good agreement with 
other estimates in literature. 

As a further test, we compute the positions of the first 
three peaks in the CMBR anisotropy spectrum using the 
procedure detailed in ^^2^. According to this prescrip- 
tion, in a flat universe made out of a matter term and a 
scalar field -like fluid, the position of the m-th peak is 
given by : 



Im = lA{m -if - Sifim) 



(24) 



with Ia the acoustic scale, (f the overall peak shift and 
Sip„i the relative shift of the m-th peak with respect to 
the first. While (p and SiPm are given by the approx- 
imated formulae in Ref. |4g, the acoustic scale for flat 
universes may be evaluated as : 



Cs 



F 



I ais 



(25) 



with : 



da 



-1/2 



(26) 



where we use Ea. (|21ll to determine a;^ = (l + zig) ^. The 
other quantities entering Eas. H25|l and H26|l are defined as 
follows Ellii : 
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FIG. 3: Constraints on the deceleration parameter go and the scaled density parameter y for the MB (left panel), Dt (central 
panel) and PR (right panel) EoS. Models with parameters on the right of the short, solid and loong dashed lines give values of 
('i, hjh) respectively in agreement within la with the measured ones. The black dot individuates the best fit model discussed 
in the text. For all the EoS, Qm is set to the best fit value reported in Table I. 
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(27) 



(28) 



(29) 



where t = / a~^dt is the conformal time, pt and pr are 
the energy densities of the baryons and radiation respec- 
tively, w{z) and fl^f, = Pcj,/ Pcrit{z) are the barotropic fac- 
tor and the density parameter of the scalar field. In order 
to use Eas. (|25|l - H29I) . we note that the role of the scalar 
field fluid is played by the dark energy so that all the 
quantities with the subscript (j) have now to be evaluated 
using the energy density corresponding to a given EoS. 
Finally, we set the present day value of the radiation den- 
sity parameter as VLq = 9.89 x 10~^ and n = 1 as 
index of the spectrum of primordial fluctuations, entering 
the approximated formulae for Cp and SLpm- 

The position of the first two peaks in the CMBR 
anisotropy spectrum ha been determined with great ac- 
curacy by WMAP giving |^ : 

;M/MAP^ 220.1 ±0.08 , = 546 ± 10 , (30) 

while the position of the third peak is more uncertain 
and may be estimated as : 



Table I for each EoS and use a analysis to constrain 
go and y. Formally, the best fit parameters turn out to 
be: 



(-0.972, 4.80) for the MB EoS , 
(90, y) = { (-0.950,0.88) for the Dt EoS , 
(-0.988, 0.34) for the PR EoS , 



(32) 



giving {li,l2,h) = (208.3,546.3,857.2) as best fit values 
independent of the EoS considered. However, as Fig.|2| 
shows, the region of the parameter space {qo,y) that 
are consistent within la with the bounds from the posi- 
tion of the peaks is quite large, so that, as already pre- 
dicted, only weak constraints can be derived. Although 
a detailed fit to the full CMBR anisotropy spectrum is 
needed, this preliminary analysis gives encouraging re- 
sults. Indeed, considering the most stringent cut (that 
on Zi), Fig.|3| shows that it is possible to find out models 
that are in agreement with both the fit to the dinien- 
sionless coordinate distances and the position of the first 
three peaks. 



VI. SCALAR FIELD POTENTIAL 

Although the agreement with the observations is a 
valid motivation for these models, it is nonetheless im- 
portant to look for a theoretical approach to further sub- 
stantiate our proposal. Such a scheme may be easily 
recovered in the framework of scalar field quintessence. 
In such a case, the energy density and the pressure of the 
dark energy fluid read : 



Boom 



851 ±31 



(31) 



Having only three data points, it is clear that only qual- 
itative constraints can be imposed on the model param- 
eters. For this reason, we fix f^M to the best fit value in 



P0 = ^</'' + ^(0) 



02 _ V{cj>) 



(33) 



(34) 
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FIG. 4: Reconstructed scalar field potential over the redshift 
range (0, 10) for the models with the MB (short dashed), Dt 
(solid) and PR (long dashed) EoS. The potential is normal- 
ized to be 1 at z = 0, while, on the abscissa, (/> is in units of 
the Planck mass Mpi. 



where is the scalar field evolving under the action of 
the self- interaction potential V{4>). For a given V{4)), 
Eas. H33|l and l34() may be inserted in the Friedmann equa- 
tions in order to determine p,f,{z), w^{z) = p^{z) / p^{z) 
and the other dynamical quantities of interest. It is worth 
noting, however, that this procedure may be reverted 
so that, for given w(z) and E(z), one may find out the 
self- interaction potential V{(f)) giving rise to that kind 
of cosmological expansion. To this end, one may first de- 
termine d(l){z)/dz and V{z) as ^3 (but see also ^3 for 
reconstruction from the SNela data directly) : 



1 



Viz)^-il~w^)Eiz) , 



dz 



1 



nM{i + z) 



2 1 



(1 ~ nM)Eiz) 



-1/2 



(35) 



(36) 



where we have defined V = V/p^{z = 0) and ~ (f)/Mpi 
with p^{0) = (1 - flAi)pcrit and Mpi = (STrG)"^/^. Note 
that, to get Ea. l(^ . we have chosen (j> > (which gives 
d(j)/dz < 0) without any loss of generality. Ea. l|5^ gives 
V{z), while V{(j)) may be obtained integrating Ea. l|5^ 
with the initial condition 0(z = 0) = to get (/'(z) and 
then inverting this relation with respect to z. 

We have applied this procedure imposing w^{z) = 
Wi{z) (with i = MB, Dt, PR) over the redshift range 
z = (0, 10) in order to recover the scalar field potential 
giving rise to our exotic EoS. Not surprisingly, an ana- 
lytical solution is not possible so that we have resorted 
to numerical techniques setting the model parameters to 
their best fit values. The results are shown in Fig.0|for 
the MB, Dt and PR EoS. It is worth noting that, al- 
though the three EoS are quite different, the potential 
V{(j)) is remarkably similar and any difference could be 
hardly detected over a large range in 0. As a matter 



fact, the same analytical approximating function may be 
fitted to the three models. Indeed, we find that, within 
2%, a very good fitting is obtained for : 



Viq^) ^VoU + Vsl^] exp 



(37) 



with (Vs, (j)s, a) fitting parameters to be determined on a 
case by case basis. For the best fit models, we get : 



{Vs,4>s,a) 



(0.0448, 
(0.0331, 
(0.0369, 



-0.1817,0.6098) 
-0.1611,0.5197) 
-0.1517,0.9434) 



(38) 



for the MB, Dt and PR EoS respectively. Considering 
the behaviour oi 'w{z) at high z shown in Fig.|21 it is 
somewhat surprising that the same functional expression 
approximates well the potential V{(j>) for all the three 
EoS. However, we have checked that there are no sys- 
tematic errors in the reconstruction procedure. Indeed, 
for the Dt and PR EoS, at high enough z, (jP is negligible 
with respect to V{(f>) so that the scalar field enters in a 
slow roll like regime and w ~ —1 as in the lower pan- 
els of Fig. 121 For the MB case, a slow roll regime is not 
achieved at high z where, on the contrary, 0^ ~ ^{4') and 
hence the EoS of the scalar field counterpart vanishes. 

It is worth noting that the approximating potential is 
quite different from those often used in literature, such as 
the exponential potential [i^ l50l | and the power law one 
[5ll |. On the other hand, its shape is the same as those 
proposed in supergravity iiispired models according to 
which it is 1/ cx (j)" exp {cjp) [53 ■ However, a < in such 
models, while we find a positive. While for large 0/0s 
both our approximating potential and SUGRA - like ones 
are exponential, for small (j)/4>s, V{(l)) takes a power law 
shape in the SUGRA case, while, for our models, V{(j)) is 
approximately constant so that a cosmological constant 
behaviour is achieved. This is consistent with the results 
w{z = 0) ~ —1 we find for the present value of the EoS 
using the best fit parameters in Table I. 



VII. CONCLUSIONS 

Notwithstanding their (somewhat radically) different 
approaches to the dark energy puzzle, all the models pro- 
posed so far assumes that the dark energy EoS is a linear 
function of the energy density. From a thermodynamical 
point of view, the ansatz p — wp (no matter whether w 
is a constant or a function of z) means that the fluid is 
modeled as a perfect gas. Putting forward the analogy 
with classical thermodynamics, however, it is well known 
that the perfect fluid approximation is quite crude and it 
is, in particular, unable to deal with critical phenomena 
such as phase transitions and the behaviour of the fluid 
near critical points. As a matter of fact, the perfect fluid 
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approximation works only in very particular conditions. 
On the other hand, our deep ignorance of the fundamen- 
tal properties of the dark energy nature does not motivate 
the choice of such idealized condition for the present state 
of this component. It is thus interesting to consider EoS 
that are more general than the perfect fluid one reducing 
to this latter in a certain regime. 

Motivated by these considerations, we have investi- 
gated here the consequences of abandoning the perfect 
fluid approximation on the dynamics of a dark energy 
dominated universe. To this end, we have considered 
four different EoS, namely the Redlich - Kwong Eq.(|SJl, 
the Modified Berthelot Eq.©, the Dieterici Eq.Q and 
the Peng - Robinson Eq.® parametrizations. These have 
been chosen because, from classical thermodynamics, we 
know they are well behaved also in critical conditions. 
The viability of the models and constraints on their char- 
acterizing parameters have been studied by using a like- 
lihood analysis taking into account the observations on 
the dimensionless coordinate distance to SNela and ra- 
diogalaxies and priors on the acoustic peak and shift pa- 
rameters A and TZ. This test has shown that all the four 
EoS are able to give rise to models that fits quite well 
the available dataset, but the RK EoS has to be rejected 
since it does not give rise to a deceleration phase in the 
past. On the other hand, the MB, Dt and PR EoS pre- 
dict reasonable values for the transition redshift zt and 
the age of the universe to in good agreement with previ- 
ous model independent estimates. As a further check, we 
have also evaluated the position of the first three peaks 
in the CMBR anisotropy spectrum finding out that, for 
each model, there exists a region of the parameter space 
such that both the CMBR peaks and the r{z) diagram 
are correctly reproduced. These successful results are 
quite encouraging since they show that the perfect fluid 
EoS may be given off without worsening the agreement 
with the data. This nice consideration strongly moti- 
vates further test of these models in order to both better 
constrain their parameters and try to select among them 
according to what model is better suited to describe what 
is observed. Since the MB, Dt and PR EoS evolve with 
z in different ways, it is desirable to resort to observ- 
ables depending on w{z) rather than its value over only 
a limited redshift range. Interesting candidates, from this 
point of view, are the full CMBR spectrum (not only its 
peaks position) and the growth factor that also depends 
on the theory of perturbations. 

As an interesting byproduct of the likelihood test, we 
have discovered a degeneracy with both the concordance 
ACDM and the quiessence (dark energy with constant 
w) QCDM models. Actually, we have fitted both ACDM 
and QCDM models to Gold -like SNela dataset simu- 
lated using one of our EoS as true background cosmolog- 
ical model. Using the same procedure in Riess et al. 
0, we have found that the ACDM model provides a 
very good fit for values of in well agreement with 
the input ones. On the contrary, the QCDM model still 
gives a good match with the simulated data, but ^Im 



is slightly biased high and w may be artificially pushed 
in the phantom region w < —1. This suggests the in- 
triguing possibility that phantom models turns out to be 
the best fit to SNela data only because of a systematic 
error on the EoS. It is worth stressing, however, that 
this result is only preliminary being based on a limited 
dataset. To further substantiate it, one should carry out 
a Fisher matrix analysis or detailed Montecarlo simula- 
tions also taking into account other probes as the full 
CMBR anisotropy spectrum. 

Having been inspired by classical thermodynamics, the 
EoS considered are phcnomenologically motivated, but 
lacks a background theoretical model. To overcome this 
difficulty, we have worked out a scalar field interpreta- 
tion reconstructing the self - interaction potential in such 
a way that the quintessence EoS is the same as the MB, 
Dt or PR EoS. The potential V{(l)) thus obtained is very 
well approximated by an analytical expression that is nei- 
ther exponential- like nor power law- like, but has for- 
mally the same form as the SUGRA inspired models. 
However, there is a significant difference since, for small 
<j>, both potentials scale as cj)", but a > for our models 
rather than negative as in the SUGRA scenario. It could 
be interesting to work out the consequences of such a dif- 
ference. On the other hand, it is also possible that the 
EoS considered have to be considered as effective ones 
as is the case in some braneworld inspired dark energy 
models or for the curvature fluid [illMIIll in f{R) 
theories. 

We would like to conclude with a general comment. 
As it is well known, the perfect fluid EoS is only a crude 
approximation of a real fluid that is usually used in cos- 
mology since it represents the simplest way to fit the 
available data. However, as first shown in the Van der 
Waals quintessence scenario |2^ [s^] and further demon- 
strated here, abandoning the perfect fluid EoS still makes 
it possible to fit the available astrophysical data with the 
same accuracy so that the use of realistic EoS turns out 
to be motivated also a posteriori. In our opinion, the era 
of precision cosmology calls for precision theory so that 
time is come to abandon approximate description such 
as the perfect fluid one. Which is the most realistic de- 
scription of the dark energy term is a topic that worths 
be addressed with the help of hints coming from thermo- 
dynamic analogies. 



APPENDIX A: SOME DETAILS ON THE EOS 

Although inspired by classical thermodynamics, the 
EoS we have considered are somewhat exotic so that we 
believe it is useful to give some further details from the 
thermodynamic point of view. As a preliminary step, 
let us denote with {p, V, T) the pressure, the volume and 
the temperature of the fluid and with a subscript c these 
quantities evaluated at the critical point. Let us also 
remember that the critical point is defined by the condi- 
tions : 
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dp 
dV 



T=cst 



fcPp_ 



= . 



(Al) 



T=cst 



Let us also denote with the subscript r the reduced quan- 
tities, i.e. Xr = x/xc- In what foUows, we show how 
Eqs.Q - (HI E^re obtained starting from their thermody- 
namical analog. As a general remark, note that the tem- 
perature of the dark energy fluid should be intended as 
an effective rather than a physical one in order to avoid 
problems with negative values. Note that this problem 
is also present in the case of the perfect fluid EoS where 
a negative u; < is formally equivalent to a negative T. 



1. Redlich - Kwong 

Let us first consider the RK EoS that is defined as : 



P 



RT 



V -b V{V + b)T^/^ 



(A2) 



with R the gas constant and (a, b) two model parameters. 
Inserting Eq. ljA2|l into Eq. IjAljl , we get : 



a=\lz- 2V2 RT^''^Vc , 6 = (1 - ^/2)Vc 
so that Eq. (|A2|I may be rewritten as : 



p = PcX 



(1 - \/2)-iv; - 1 I V3-2V2 



1 



(A3) 



with : 



2(3 - 2V2) 



(l-y2)(l- V3-2V2) 



X Pc , 



V 2 Vc 

Note that the critical pressure is almost half the perfect 
fluid one. Let us now assume that the temperature is 
constant and equal to its critical value^ so that Tr = 1. 
Using then V = 1/p ^ Vr = Pc/ P, after some algebra, 
Eq. (|A3p may be finally rewritten as : 



p — Pt] X 



1 - \/3 ~ 2V2ar] 
1 - (1 - V2)ari 



^ The ansatz T = Tc is somewhat arbitrary, but does not lead 
to any loss of generality. Actually, choosing a different T only 
rescales the EoS. Since we do not know either T or Tc, it is a 
useful working hypothesis to set Tr = 1. 



which is the same as Eq.® having posed p = p/pcrit, 
rj = p/ Pcrit and defined : 



a = Pcrit/ Pc 



^/2 



1 - ^3 - 2V2 



X Pc/Pc 



Note that, for the best fit parameters in Tablet, /3 < 
so that < as expected since dark energy is known to 
have negative pressure. 



2. Modified Berthelot 

Let us summarize here the main steps above for the 
case of the MB EoS. The pressure p as function of tem- 
perature T and volume V is implicitly defined as : 



p : 



RT 



1 



p 



128 \p. 

where the critical pressure is given as : 

83 RTr 



1 



y2 

^ C 

2^2 



(A4) 



Pc 



128 Vc 

Using this relation, we may rewrite Eq. ljA4|l as : 



128 

P = Pr X 

^ 83 ^ 



Tr/Vr 



1- (9/128) V;r^(l-6r,r^) 



(A5) 



Setting Tr — I and Vr — Pc/ P, after some simple algebra, 
we finally get : 



P 



1 arj 

with p — p/ Pcrit 1 f] = p/ Pcrit and we have defined : 



128pc 



/3=t^^-4.4(WPc). 
45 X 83 Pc 



As for the RK case, the best fit parameters in Table 1 
gives /3 < as expected. 



3. Dieterici 

In terms of thermodynamical quantities, the Di EoS 
reads : 



RT 

V-b 



exp 



RTV 



(A6) 



where the two parameters (a, b) may be determined solv- 
ing Eq. (|Aip thus obtaining : 



o = 2RTcVc 



C' C •) 



b = Vc/2 . 
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Introducing reduced variables {Tr,Vr) and using the ex- 
pression for the critical pressure : 

2i?T, 



we rewrite Eq. (|A6p as : 



P 



PcTr 

2Vr - 



-exp[2(l-T-iv;-i] . (A7) 



With the usual positions = 1 and Vr = PcP, the above 
relation finally becomes : 



P = 



f3v 



2 — arj 



exp [2(1 — arj)] 



which is the same as Eq.Q provided one defines tilted 
quantities as usual and 

a = Pcrit/P , /3 = Pc/Pc ■ 

Not surprisingly, /3 (and hence Pc) turns out to be nega- 
tive for the best fit parameters. 



4. Peng - Robinson 

As a final case, let us consider the PR EoS starting 
from its expression in terms of thermodynamical quanti- 
ties : 



RT 



V-b V{V + b) + biV -b) 



(A8) 



Solving the equations for the critical points allows us to 
express (a, b) in terms of (Tc, Vc) giving : 



a = CaRTcVc , b = CbVc 

with Ca — 1.487 and Cf, ~ 0.253. The critical pressure 
turns out to be : 



RTc/Vc 

1 - Cf, 



{l+Cb)/{l-Ci,)+Cb 



0M7RTjVc 



Introducing reduced variables, Ea. (|A8p rewrites as 



PcTrCp ^ 
Vr - Cb 



1 - 



Vr{Vr+Cb)l{Vr~Cb)+Cb 



(A9) 



Eq. (|HJ| is finally obtained by setting in the above relation 
Tr ~ 1 and Vr — Pc/p thus getting: 



P = 



1 — arj 



1 



{ca/cb)ari 



(1 + ar])/{l — ar]) + arj 
with p and rj defined as usual, while we have set : 

Pc/Pc 



OL^PcritjPc , P 



CbCp 



Once again, the best fit parameters gives /3 < and hence 
Pc < as expected. 
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